Solitons in a chain of PT-invariant dimers 
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Dynamics of a chain of interacting parity-time invariant nonlinear dimers is investigated. A dimer 
is built as a pair of coupled elements with equal gain and loss. A relation between stationary soliton 
solutions of the model and solitons of the discrete nonlinear Schrodinger (DNLS) equation is demon- 
strated. Approximate solutions for solitons whose width is large in comparison to the lattice spacing 
are derived, using a continuum counterpart of the discrete equations. These solitons are mobile, 
featuring nearly elastic collisions. Stationary solutions for narrow solitons, which are immobile due 
to the pinning by the effective Peierls-Nabarro potential, are constructed numerically, starting from 
the anti-continuum limit. The solitons with the amplitude exceeding a certain critical value suffer 
an instability leading to blowup, which is a specific feature of the nonlinear PT-symmetric chain, 
making it dynamically different from DNLS lattices. A qualitative explanation of this feature is 
proposed. The instability threshold drops with the increase of the gain-loss coefficient, but it does 
not depend on the lattice coupling constant, nor on the soliton's velocity. 

PACS numbers: 42.25.Bs, ll.30.Er, 42.82.Et, 42.81. Qb 



I. INTRODUCTION 

In the original works by C. Bender and co-authors 
[l|, , it was pointed out that non-Hermitian Hamiltoni- 
ans can have an entirely real eigenvalue spectrum under 
the parity-time (VT) symmetry constraint. This mathe- 
matical observation can have deep physical consequences, 
essentially altering the familiar properties of collective 
modes in the respective media combining amplification 
and dissipation. Indeed, it is usually assumed that the 
balance between the gain and loss uniquely selects pa- 
rameters of isolated stable modes. However, the use of 
the VT symmetry makes it possible to support continu- 
ous families of modes, allowing the dissipative media to 
emulate conservative ones, up to a certain threshold level. 
In particular, the effects of amplification and dissipation 
may stay in the exact balance for symmetric configura- 
tions of weak fields, hence the 'PT'-invariant dynamics 
is preserved. In contrast, for field intensities above the 
threshold, nonlinear self-action (if present in the system) 
breaks the VT symmetry both locally and globally, re- 
sulting in the asymmetric wave localization in the region 
with amplification. 

As the necessary condition for the T'T-symmetry of 
the Hamiltonian with a complex potential, V{x), is re- 
duced to condition V{x) = V*{—x), such systems can 
be realized in the most straightforward way in optics, by 
combining a spatially symmetric profile of the refractive 
index with symmetrically placed mutually balanced gain 
and loss [j] (complex potentials may also be physically 
relevant in the case when they are not subject to the 
T'T-symmetry constraint Q). The possibility to realize 
physical systems with the VT symmetry was a motiva- 
tion for many theoretical [H and experimental works. 



Effects of the conservative nonlinearity in PT-symmetric 
systems were addressed too 043 • Recently, the nonlin- 
earity of the gain and loss, also subject to the condition 
of the VT symmetry, was introduced in Ref. ■ 

The simplest nonlinear T'T-symmetric object can be 
realized as a pair of linearly-coupled optical waveguides 
(alias a dimer), composed of a passive waveguide carry- 
ing linear loss and its active counterpart imparted with a 
matched compensating gain |8] (in a more general form, 
without the condition of the exact equilibrium between 
the loss and gain, the same system of linearly coupled ac- 
tive and passive waveguiding cores was investigated for 
a long time in various forms [iH-[i3, see also a review 
in Ref. [16| ) . A chain of such PT- invariant couplers was 
studied in Ref. @, with each active or passive element 
linearly connected with an element of the opposite sign, 
belonging to an adjacent coupler. In other words, in this 
setting the axis of each waveguiding coupler is perpendic- 
ular to the direction along which the chain is built. It 
was demonstrated that this chain of couplers can support 
stable solitons with amplitudes smaller than a threshold 
amplitude. 

In the present work, we propose and analyze another 
array of T'T-symmetric couplers with the intrinsic non- 
linearity. In the array, active and passive elements are 
linearly coupled to the elements of the same sign belong- 
ing to adjacent dimers [see Fig. [Ha)], i.e., the axes of 
the couplers are aligned with the direction of the chain 
composed of them. In fact, the proposed system is quite 
natural in terms of the realization in optics, as it can be 
built of two parallel extended waveguides, one pumped 
and one lossy, if each one is segmented into an array of 
individual waveguides. Our aim is to find nonlinear lo- 
calized modes (discrete solitons) in this chain, and study 
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their stability, mobility, and interaction. 

The paper is organized as follows. Section ITTl outlines 
the model. Analytical results for moving and stationary 
solitons and their stability are presented in Sec. IIIII Re- 
sults of numerical studies are reported in Sec. IIVI and 
Sec. IVl concludes the paper. 



II. MODEL 

As said above, we consider a nonlinear chain shown 
in Fig. [T{a), which is described by the following set of 
equations: 

at 

iC (m„+1 + Un-l - 2Un) , 

zC(z;„+i 2w„) . (1) 

Here 7 > is the coefficient of the gain and loss act- 
ing on complex variables u„ and u„, which correspond, 
respectively, to the active and passive elements, real co- 
efficient a accounts for the nonlinear frequency shift (the 
intrinsic nonlinearity of each element), n is the constant 
of the vertical coupling inside the dimer, and C is a coef- 
ficient of the horizontal linear coupling between dimers in 
the chain, see Fig. [Ija). Using obvious transformations 
{vn — > —Vn, the complex conjugation, and staggering, 
respectively), one can fix k, C, and a to be positive. Fi- 
nally, by means of rescaling, we can set k = cr = 1, which 
leaves 7 > and C > as irreducible parameters, the 
final form of the equations being 

= ^Un+i\Un\'^Un+iVn+iC{Un+l+Un-l-2Un)., 
= --iVn + i\Vn\'^Vn+iUn + iC{Vn+l+Vn-l-'2Vn). 

(2) 
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Figure 1: (a) A schematic of the chain of dimers. (b) The 
spatio-temporal evolution of |u„(t)p showing the propagation 
of a soliton through the chain of dimers described by Eq. ((2]) 
[the respective plot for is nearly identical to this one]. 

The initial conditions were set as per the approximate solution 
dSl), with 7 = 0.5, A = 0.37, C = 0.5, V = 0.2. (c) Same as 
in (b), but for a larger velocity, V — 2. 



First, we will derive approximate solutions for mov- 
ing broad solitons, whose width is large in comparison 
with the lattice spacing, and study their almost elas- 
tic collisions, using a continuum counterpart of the dis- 
crete equations. We will also consider the case of narrow 
solitons, finding such strongly localized solutions numer- 
ically, starting from the ant i- continuum limit. It will be 
demonstrated that the solitons with the amplitude ex- 
ceeding a certain critical value become unstable. Because 
the instability is induced by the interplay of the gain and 
loss in the T^T-symmetric system, the instability thresh- 
old naturally lowers with the increase of the strength of 
the gain and loss. 



III. SOLITONS 

A. Analytical approximations 

We start with the analysis of the continuum limit (C — > 
00) of Eq. (12), which leads to the system of PDEs, 

du I i2 • 

_ = +^u + l\u\ u + + 

(^"^ -I l2 • d^"" /ON 

— = -^v + i\v\ v + iu + i—^, (3) 
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which possess exact soliton solutions, 
A 



u{t, x) 
v{t, x) 



A sech 



Here A, A ^, w, and 5 represent the ampUtude, width, 
velocity, frequency, and intrinsic phase shift of the soli- 
ton: 



~A^ - cos J + V'^/A + ^2/2, sin 8 



-7, 



cos(5 ±r,r= 



r 



(5) 



In the following the case cos (5 = F will be considered. 
Note that T < 1, and this solution exists for 7 < 1. In 
the limit of the vanishing gain and loss, 7 — > 1, Eqs. ([3]) 
become tantamount to the well-known model of dual- 
core nonlinear optical fibers. In that case, soliton ^ 
with cos (5 = r = +1 reduces to the symmetric soliton 
in the dual-core fiber, which loses its stability through 
the symmetry-breaking bifurcation at a finite value of 
the energy [13]. For the comparison with the discrete 
system introduced in the present work, it is important to 
mention that a similar symmetry-breaking instability of 
two-component solitons in the discrete counterpart of the 
dual-core- fiber model (in fact, this discrete system may 
also be readily implemented in optics) was found in Ref. 

Although equations dU are dissipative, they are 
Galilean invariant, hence the soliton may move at an ar- 
bitrary velocity, V . It is also worthy to note that the am- 
plitude of the soliton may be arbitrary, i.e., the present 
solutions form a continuous family, unlike formally sim- 
ilar approximate [ll|, [3 Ell and exact [l^ soliton solu- 
tions previously found in the above-mentioned models of 
dual-core systems with the gain applied in one core, and 
loss acting in the other. Those models were considered 
without the condition of the balance between the gain 
and loss in the general case (although the case of equal 
gain and loss had a special purport in that context too), 
therefore they support solitons with two discrete values 
of the amplitude, one stable and one unstable. 

For sufficiently wide solitons, i.e., for A <^ VC, the 
following approximate solution to the discrete equation 
can be deduced from the continuum solution @, 



Unit) 
Vn{t) 

X exp 



A sech 



A 
71 



Xq 



A^ y2 



c 



t + 



vt 



-iS/2 



V{n - xa) 



2VC 



(6) 



where xq is shift along the lattice. An example of the 
soliton propagation, generated by simulations of Eq. ^ 
with the initial conditions taken as per approximation 
([6|), is displayed in Fig. [Ijb). In this case, A/\/2C — 
0.37, i.e., the soliton is not very wide; nevertheless, ap- 
proximation © yields good initial conditions for the soli- 
ton. The moving soliton practically does not radiate en- 
ergy in the form of small-amplitude waves. On the con- 
trary. Fig. [IJc) demonstrates that the soliton with a large 
velocity lose energy through the emission of radiation. 



For narrow solitons, with A > v C, the approximation 
based on the continuum limit is obviously irrelevant. In 
this case, solitons can be constructed numerically, start- 
ing from solutions which are exact ones in the anticon- 
tinuum (AC) limit, that corresponds to C = in Eq. 
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= A 



gi<5/2 

for 



(7) 



where A is an arbitrary amplitude, while the frequency 
and phase shift between the two components are given 
by 

uj — -A^ ~ C0S.5, sin (5 = —7, cos(5 = ±r. (8) 

In the following the case cos (5 = F will be considered, 
following the choice adopted above for the broad solitons. 
It follows from Eq. ((8)) that this solution does not exist 

at 



7 > 1 or u < — F. 



(9) 



On the other hand, for C ^ 0, one has the continuous- 
wave (CW) solution to Eq. ([2]), i.e., the wave with a 
constant amplitude, in the form of 

W„ = UO,Wn=WO (10) 

for all n, where uq and vq are defined by Eq. ((T]) and Eq. 
®. 

We stress that, as well as the exact soliton solutions 
(j4|) and approximate ones ([6]), the AC and CW solutions, 
given by Eqs. ([7]) and pH)) . respectively, form continuous 
families in spite of the fact that they exist in the dissipa- 
tive system. This is a fundamental manifestation of the 
VT symmetry of the present system. 

B. Stability analysis 

The next step of the analysis is the consideration of the 
stability of the AC solution (O and CW solution given 
by Eq. pil| . In the latter case, the perturbed CW state 
is looked for as 



r {A+ 



eie 



at — ikn 



at~ikn 



(11) 



where ei and £2 are complex amplitudes of infinitesimal 
perturbations with eigenvalue a and wavenumber k. Sub- 
stituting this expression into Eq. ([2]) and performing the 
linearization with respect to the perturbations yields four 
branches of the dispersion relation: 



ai,2 = ±2y/C{l - cos k) [v42 - C(l cos k)] 

= ±2VX, (12) 
asA = ±2^/x + r [-r + ^2 _ 2C(1 - COS fc)]. (13) 
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The CW solution is stable if condition Re{a} < holds 
for all four eigenvalues at all k. For C ^ it is suf- 
ficient to analyze this condition for 0:3.4, as inequality 
Re{Q;i,2} < follows from Re{Q;34} < 0. The subse- 
quent consideration demonstrates that the CW solution 
might be stable only if condition 



< C(l -cos/c) 



(14) 



holds for all k, which is impossible at A 0. Thus, 
for C 7^ 0, all the CW solutions pU)) are modulationally 
unstable. This result also suggests that dark solitons are 
unstable in the dimer chain, because the corresponding 
CW background cannot be stable. 

However, in the case of periodic boundary conditions, 
rather than the infinite chain, the CW state can be stable 
for C 7^ if the ring-shaped chain is sufficiently short. 
As follows from Eq. (ITi)) . unstable are the modulational 
perturbations with wavelengths 



A 



2tt 



2n 



> 



arccos(l - A'^/C) ' 



Thus, if the number of dimers in the chain with periodic 
boundary conditions satisfies condition 



N < 



2-K 



arccos(l- A2/C)' 



(15) 



which makes the length of the chain smaller than the 
wavelength of the shortest unstable perturbation,the CW 
solution ([TU| is stable. Because N cannot be smaller than 
2, a constraint 



A^ < 2C 



(16) 



follows from Eq. (jTS]). 

The stability analysis for the AC solution Eq. ([7]) can 
be performed by substituting Eq. (|lip into Eq. ([2]) with 
C = 0. In the case, Eqs. (fT2|) and (fT3l) reduce to ai^2 = 



and a3_4 — ±2-\/r (— F -|- A^), from where it follows that 
the AC solution is stable for 



A. Unstable CW dynamics in the chain with 
periodic boundary conditions 

Interesting dynamics was observed in the chain with 
periodic boundary conditions for the initial conditions 
taken as per the CW solution (fTO)) . in the case when 
the chain's length, N, is only slightly greater than the 
critical value specified by Eq. (IT5|) , while conditions 
and pri) are satisfied. 

The instability of the system results in periodic appear- 
ance and disappearance of weakly localized states. This 
dynamics is illustrated by Fig. [5] for 7 = 0.5, A = 0.3, 
C — 0.3, N — 8 (the corresponding CW solution is sta- 
ble for N < 7). Panel (c) in Fig. [5] displays the time 
variation of the localization parameter defined as 



I'JV-l 



El 

\n=0 



(18) 



Note that when = const, as in the unperturbed CW 
state, Eq. (|18l) yields / = 1/N, which is the minimal 
possible value, which tends to zero with the increase of N. 
On the other hand, when single |u„| is different from zero, 
the localization parameter attains its maximal value, I — 
1. 

In these simulations, the only source of perturbations 
destabilizing the CW state (ITU)) was rounding errors of 
the numerical code. In this case, the deviation of the 
solution from the initial unstable configuration becomes 
visible at t « 900. The unstable dynamics a.t t > 900 
results in quasi-periodic localizations and delocalizations 
of |u„p, as seen in Fig. [2jc). The maximum of |u„p 
gradually moves along the chain with a nearly constant 
velocity, as seen in panels (a) and (b). 

A similar effect of the periodic in time and space local- 
ization of the unstable CW solution, subject to periodic 
boundary conditions, was observed in the model of a di- 
atomic chain of atoms To']. A simple analytical model 
describing this effects was proposed in Ref. [20j . 



B. Strongly localized on-site solitons 



A'< Ai?2 



(17) 



These solutions, stable in the AC limit, are used below 
to construct numerical solutions for narrow solitons. 



IV. NUMERICAL RESULTS 

Discrete evolution equations ^ were integrated nu- 
merically, using a modification of an implicit, uncondi- 
tionally stable Crank-Nicholson scheme with accuracy 
~ r*, where the time step was taken as r = 2 x 10~^. 
Zero boundary conditions were employed except for Sec. 
IIV Al where the periodic boundary conditions were used. 



Sharp (tightly localized) on-site solitons were con- 
structed starting from the AC- limit solution (O, which 
is valid for C — 0. After setting the initial conditions, 
the coupling constant was increased linearly in time from 
C = 0att = 0toC== 0.15 at t = 400, and then kept 
constant. The total duration of the numerical run was 
t = 1000. With such a slow increase of C, nearly station- 
ary solitons were readily created (details of the station- 
ary shape of the solitons are given below). The initial 
amplitude of the soliton in the AC limit was A — 0.6. 
This amplitude satisfies the AC stability condition ([T7]) 
for three considered values of the gain/loss parameter, 
7 = 0, 0.3, 0.9. The chain was composed of 21 dimers, 
which was sufficient for constructing tightly localized soli- 
tons. 
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Figure 2: Unstable dynamics in a compact (A'' = 8) ring- 
shaped chain of dimers with initial conditions corresponding 
to the CW state (|10[) with amplitude A = 0.3 and model's 
parameters 7 = 0.5, C = 0.3. (a) Sites (dimers) in the (n, t) 
plane with intensity greater than I.IA^ are shown in black, 
while others are not shown, (b) The evolution of intensity 
[the picture for is identical to this one, within 

the numerical accuracy] . (c) Evolution of the localization pa- 
rameter, I, defined as per Eq. (|18[) . 



Results of the simulations can be summarized as fol- 
lows: 

(i) The loss-gain coefficient 7 does not affect the pro- 
file and amplitude of the soliton, |u„P, |wnP- This can 
be seen in Fig. [31 where the shapes of the solitons con- 
structed for 7 = 0, 0.3, 0.9 practically overlap. Coeffi- 
cient 7 only affects the phase difference between complex 
quantities u„ and u„ and the frequency of the soliton. 

The fact that the profiles of stationary solitons do not 
depend on 7 can be explained by noting that if y„ = 



is a stationary solution of the standard discrete 



0.14 



0.07 



0.00 




Figure 3: Soliton profiles obtained from the anti-continuum 
limit by the slow increase of coupling constant C at three 
different values of the gain-loss parameter, 7 = 0, 0.3, 0.9. 



Intensities \u„\ and 



do not depend on 7, while the soli- 



ton's frequency and phase shift between Un and Vn depend on 
7, as explained in the text. 



then a stationary solution of Eq. ([2]) can be written as 



pf5/2 
~iS/2 



(20) 



where 6 is defined by Eq. It is clear that the solu- 
tion (l^ni) with 7 7^ differs from a stationary solution 
to Eq. (|I9I) . with 7 = 0, only by the frequency and 
phase shift, while their intensities are same. Thus, the 
discrete system based on Eq. ([2]), although being dissipa- 
tive, supports the continuous family of solitons, being in 
that respect similar to the conservative DNLS equation, 
and different from generic discrete dissipative systems, 
such as the discrete complex Ginzburg-Landau equation 

(ii) To check that the numerically found solitons are 
stationary ones, we calculated the measure of non- 



stationarity at site n, S„ 



/dt for three soli- 



nonlinear Schrodinger (DNLS) equation, 



tons constructed from the AC limit at 7 — 0, 0.3, 0.9. 
We found that the result does not depend on 7, as ex- 
pected in view of the above comment. In the considered 
examples, the non-stationarity of the constructed soli- 
tons is weak (max{5'„} does not exceed 0.0003), and it 
gets weaker for the solitons constructed with the use of 
C increasing still slower in time. 

(iii) The tightly pinned on-site-centered soliton re- 
mains stable if its amplitude is not too large to break 
the VT symmetry. This conclusion is justified by mon- 
itoring the long-term dynamics of the discrete solitons 
constructed according to the protocol outlined above. 



Instability of inter-site and twisted tightly 
localized solitons 



dyn 
dt 



-iyn + i\yn\^yn + iC{yn+l + J/n-l - 2?/„), (19) 



On the contrary to the stable on-site centered discrete 
solitons, outlined above, the analysis has revealed that 
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Figure 4: The evolution of the inter-site-centered soliton, 
grown from the AC initial conditions, following the increase 
of C at two different values of 7. Eventually, one site spon- 
taneously sucks the energy of the other one. The parameters 
are A = 0.6, C = 0.15, N = 16. In panel (a), 7 = 0.2, and in 
panel (b), 7 = 0.6. 



all the inter-site-centered solitons are unstable. We at- 
tempted to construct such solitons by starting, in the AC 
limit, with two adjacent excited sites, n = 0, 1, using Eq. 
([7]) with A — 0.6 as initial conditions, and gradually in- 
creasing C from C = at t = to C 0.15 at i = 400, 
and keeping then C constant. 

In Fig. 0] we present the evolution of the inter-site 
mode in the course of the increase of C. This mode is 
stable only in the absence of the coupling, C — Q. At 
C 7^ 0, one site spontaneously sucks the energy from the 
other, which end us with the establishment of a stable 
on-site-centered soliton. Note that in the examples pre- 
sented in Fig. m this happens at i « 50, i.e., before C 
reaches the target constant value of 0.15. We thus con- 
clude that the inter-site soliton is unstable. This result 
seems very plausible in view of the above-mentioned re- 
lation of stationary solitons in this system to those in the 
conventional discrete NLS equation, Eq. (|19l) . where all 
the inter-site-centered solitons are unstable too [22| . 

We also attempted to excite twisted (antisymmetric) 
solitons, by starting, in the AC limit, with two adjacent 
excited sites, at n = and 1, using Eq. ([7]) with A — 0.6 
and introducing the phase shift of tt between the sites. 
This mode too turns out to be unstable at C > 0. In the 
usual lattices of the DNLS type, twisted solitons usually 
exist at large amplitudes 22], but the soliton's ampli- 
tude in the present T'T-symmetric chain is limited by 
condition ([TT]), which is a possible reason of the failure 
in looking for stable twisted solitons. 

D. Instability due to the "PT-symmetry breaking 

According to the stability condition for the soliton in 
the AC limit, given by Eq. p7)) . increase in 7 results in 

a decrease of the critical value of the amplitude, Amax, 
above which the localized mode becomes unstable. For 
the solitons in the chain with C > 0, the corresponding 
critical value was determined numerically. To this end, 
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Figure 5: Intensities of the fields at the central dimer, n = 0, 
as functions of time, for the soliton in the chain of dimers 
with a small gain-loss mismatch, e = 10~^, and 7 = 0.4, 
C = 0.15. Thick and thin lines pertain, respectively, to the 
elements with gain and loss, i.e., |«np and |unp. The abrupt 
growth of |unp signalizes the onset of the instability. 

we took, as the initial condition, a soliton constructed 
as described above (using the slow increase of coupling 
constant C), and then switched on a weak net pump 
of the energy into the soliton, by introducing a small 
mismatch, e, in the gain/loss balance. In this case, Eq. 
([2]) assumes the following form: 

I—— = +(7 + e)?i,i + I u„ Un + 
at 

iC'{Un+l + Un~l - 2Un) + W„ , 

= -{1 - ejVn+l\Vn\ Vn + 

iC{v„+l + U„-l - 2Vn) + iUn- (21) 

In the stable regime, the field intensities in all dimers 
remain equal for the active and passive elements, both 
slowly increasing due to small excess gain, e > 0. The 
system becomes unstable when the soliton's amplitude 
attains some critical value, Amax. An example of the 
unstable behavior of the central dimer (n = 0) is shown in 
Fig. [5l In the unstable regime, the intensity in the active 
element sharply increases, while in its passive counterpart 
the intensity decays. 

In Fig. |6] we compare the dependence of A^^^ on 
7, found numerically in this way, to the analytical re- 
sult for the local mode in the AC limit, (Amix)^, as de- 
fined by Eq. pT)) . It is seen that they are in a good 
agreement, despite the fact that A^^^ pertains to the 
chain with C ^ 0. The numerical value is somewhat 
higher than the analytical estimate, which may be ex- 
plained by a particular criterion of instability adopted in 
the numerical study: The soliton was assumed to be- 
come unstable when the largest (over n) relative dif- 
ference between |u„p and |wnP attained the level of 
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Figure 6: The critical amplitude, A^^^, above which the soli- 
ton's instability sets in, as a function of the loss-gain strength, 
7. The thick line is the numerical result for the solitons at 
C = 0.15. The thin line represents the critical amplitude, 

^j4max j = r, as predicted by the stability condition p7|l for 

the soliton in the anti-continuum limit (C = 0). 




max„ _ |„^|2|/ (j^^p ^ |^^j2)| > q 

Results presented in Figs. [5] and Fig. [6] were obtained 
for the inter-dimer coupling constant C ~ 0.15. Similar 
results were obtained for other values, up to C = 0.5. In 
all the cases considered, the numerically found instabil- 
ity threshold, A^axi '^^^ ^^^Y close to (Amax)^ predicted 
by Eq. p7|) for the AC-limit soliton. Thus, the thresh- 
old practically does not depend on the coupling constant 
C. This result demonstrates a drastic difference in the 
dynamical properties of the discrete solitons in the VT- 
symmetric chain from its counterparts in the usual DNLS 
equation, where the instability of various soliton modes 
is solely determined by C p^ . 

The dynamical blowup of unstable solitons is displayed 
in Fig. [71 Initial conditions were set using approxi- 
mation ([6]) with xo — 0. In this case, — 0.81 > 
r = yj\ - 72 w 0.436, hence the AC stability con- 
dition Eq. (|17p is violated. The result displayed in 
panel ^c) suggests that, for t < 40, the difference 
A = ||uo(i)P — I'yoCOPI increases with time exponen- 
tially, and then it starts to grow still faster. This change 
in the unstable dynamics is due to a qualitative change 
in the soliton's profile, that can be seen in panels [7lja,b). 
At i < 40, intensities |moP and \vo\'^ grow with time, but, 
after the difference between them becomes sufficiently 
large, |uop starts to grow faster, while |woP begins to 
decrease. 

Results presented in Fig. [7]were obtained for the quies- 
cent soliton [V — 0). We have also studied the influence 
of the soliton's velocity, < ^ < 0.5, on the onset and 
development of the instability. It was found that the 
growth rate of A in the unstable regime and the critical 
soliton amplitude A^^^ do not depend on V. 



Figure 7: The soliton dynamics in the unstable regime. Pan- 
els (a) and (b) display |ii„(t)p and |i'n(i)P, respectively. 
Panel (c) shows the difference in the intensities between the 
active and passive elements at the central dimer (n — 0), 
A — 1 1 Mop — I'^oPl, as a function of time. Parameters are 
7 = 0.9, C = 2.3, A = 0.9, V = 0. 



This instability can be explained by comparison with 
the system of two linearly-coupled DNLS equations, 
where the symmetric discrete solitons is destabilized by 
the symmetry-breaking bifurcation 18[ . It seems plausi- 
ble that the spontaneous trend to the symmetry break- 
ing, due to the self-attraction in each chain, explains 
the onset of the instability in the present setting. How- 
ever, a drastic difference from the dual-core DNLS sys- 
tem, where the symmetry breaking replaces the original 
destabilized soliton by stable asymmetric ones, in the VT 
system this is impossible, as asymmetric modes cannot 
maintain the balance between the gain and loss, and in- 
stead exhibit the blowup, as shown in Fig. [T] 



E. The Peierls-Nabarro potential 

The effect of the Peierls-Nabarro (PN) potential in the 
chain of dimers is demonstrated by simulations presented 
in Fig. [51 where the initial conditions were set by using 
Eq. 1^ with xo = 0. The parameters were chosen so as to 
produce a relatively sharp soliton when the PN barrier is 
noticeable. The soliton was kicked with a relatively small 
velocity {V = 0.02). The momentum given to the soliton 
was not sufficient to overcome the PN barrier, therefore. 



8 




Figure 8: The effect of tlie Peieris-Nabarro (PN) potential on 
tlie soliton dynamics. Tlie initial momentum imparted to the 
soliton is not sufficient to overcome the PN potential barrier, 
therefore the soliton oscillates near a local minimum of the 
potential. The parameters are A — 0.6, V = 0.02, 7 = 0.5, 
C = 0.2. 

instead of propagating along the chain, the soliton oscil- 
lates in a local well of the potential. 



F. Collisions between moving solitons 

The existence of moving solitons suggests a possibil- 
ity to study collisions between them. Examples of colli- 
sions between identical solitons with opposite velocities 
{V = ±0.3) are presented in Fig. [HI In this case, the soli- 
tons are relatively wide, hence the approximate solution 
(O can be used for setting the initial conditions. The 
collisions are seen to be nearly elastic, with the solitons 
almost completely recovering their original shapes. 

The solitons with phase shift tt in (a) repel each other, 
while the in-phase ones in (c) interact attractively. The 
collisions of the repelling out-of-phase solitons does not 
result in an overlap of their cores. For the phase shifts 
different from tt, the overlap takes place, hence the largest 
intensity at the collision point is greater than in isolated 
solitons [see panel (b)] , attaining the maximum for the in- 
phase collisions, as can be seen in (c). It is worthy to note 
that, although the instantaneous value of the amplitude 
at the collision point may be considerably larger than 
the critical amplitude defined by Eq. PT)) . this does not 
result in the onset of the instability, as the collision time 
is insufficient for that. 



V. CONCLUSIONS 

We have introduced a novel model describing an ar- 
ray of harmonically coupled 'PT-symmetric dimers. Each 
dimer consists of linearly coupled active and passive ele- 
ments with balanced gain and loss. The nonlinearity is 
represented by the conservative cubic term in the equa- 




Figure 9: Examples of quasi-elastic collisions between solitons 
for 7 = 0.5, C = 2, A = 0.5, and velocities V = ±0.3. The 
phase shift between the solitons in panels (a), (b), and (c) is, 
respectively, tt, 7r/2, and 0. 



tion of motion for each element. 

Using numerical simulations and analytical approxi- 
mations, it was demonstrated that, in the regime of the 
sufficiently strong inter-chain coupling, the weakly dis- 
crete chain supports the free motion of solitary waves. 
At sufficiently small velocities, the moving solitary waves 
do not radiate energy and collide practically elastically, 
while solitons moving with large velocities emit radiation. 

The chain of dimers was demonstrated to support 
tightly localized (sharp) stationary solitons in the regime 
of weak inter-dimer coupling (strong discreteness). Nat- 
urally, these strongly pinned solitons are immobile. A 
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relation between the stationary solutions in the chain of 
dimers and the solutions to the standard discrete nonlin- 
ear Schrodinger equation was established. 

Both wide and sharp solitons develop an instability if 
their amplitudes exceed the critical value, which leads 
to the blowup. This is a specific feature of the VT- 
symmetric chain, making the dynamical properties of the 
discrete solitons different from those in the usual DNLS 
system, in spite of the similarity in their shapes. An 
explanation to this feature, based on the trend to the 
spontaneous symmetry breaking, induced by the intrin- 
sic self-attraction in the parallel chains, was proposed. 
The critical value decreases with the increase of the gain- 
loss parameter. Above the instability threshold, the VT- 
symmetry breaking occurs, resulting in the divergence 
of the intensity at the active element, and decay of the 
intensity at the lossy one. A noteworthy finding is the 



independence of the critical amplitude on the lattice cou- 
pling constant and soliton's velocity. 

To continue the work in this direction, it would be in- 
teresting to study the chains consisting of more complex 
■PT-symmetric elements (such as those with the nonlin- 
earity of the gain and loss and to extend the anal- 

ysis for two-dimensional lattices. 
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